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We calculate parton and generalized parton distributions in Minkowski space using a scalar prop- 
agator with a pair of complex conjugate poles. Correct spectral and support properties are obtained 
only after careful analytic continuation from Euclidean space. Alternately the quark distribution 
function can be calculated from modified cutting rules, which put the intermediate state on its 
complex mass shells. Distribution functions agree with those resulting from the model's Euclidean 
space double distribution which we calculate via non-diagonal matrix elements of twist-two opera- 
tors. Thus one can use a wide class of analytic parameterizations of the quark propagator to connect 
Euclidean space Green functions to light-cone dominated amplitudes. 
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I. INTRODUCTION 



Understanding strong QCD aspects of hadron properties remains a challenge. Experimental probes at large mo- 
mentum transfer will continue to yield a wealth of data on hadron structure — stimulating theoretical explanations of 
the underlying physics in terms of quark and gluon degrees of freedom. It is well known that light-cone correlation 
functions are relevant for describing hard processes Q, since struck constituents rebound at speeds near that of light. 
In recent years there has been renewed interest in the connection between inclusive and exclusive reactions at large 
momentum transfer 0. The underlying connection is encompassed by generalized parton distributions (GPDs) which 
are functions that enter in the description of a variety hard exclusive processes 0] ■ 

Attempting to describe the non-perturbative light-cone correlations that enter in large momentum transfer processes 
has led to the formulation of gauge theories on the light cone Q. There has been progress in directly solving for 
pion light-cone Fock components in light-cone Hamiltonian QCD Another important stride has been made in 
enumerating and classifying hadronic light-cone Fock-space amplitudes gg. These could be used for modeling hadrons, 
although Lorentz covariance requires infinitely many Fock components and the incorporation of symmetries into model 
wave functions has been very limited. 

In a different approach, QCD models based on solutions to Dyson-Schwinger equations provide a useful framework 
for exploring strongly interacting bound states, see e.g. 0. This framework is fully Poincare covariant, allows for 
close contact with lattice simulations and provides a means to preserve symmetries and implement quark and gluon 
confinement. Dyson-Schwinger models have also been used to study light-cone dominated amplitudes. Calculation of 
quark distributions in the impulse approximation was undertaken in |l2f. where the nucleon Bethe-Salpeter equation 
was solved in a diquark spectator model. That investigation relied upon the use of free particle propagators for the 
quarks and diquark. In a different study, the authors of Ref. |13| calculated quark distribution functions for the pion 
using a Dyson-Schwinger type model based on entire functions. Their analysis avoided two problems: the integral 
over the relative light-cone energy is not convergent in the complex plane because non-constant entire functions are 
unbounded; secondly, expressions were derived supposing the existence of a Kallen-Lehmann representation which 
admittedly does not exist for their model propagator. Although such a model is successful at describing low mo- 
mentum, space-like processes where the propagator can be approximated as entire, the non-analytic points in the 
whole complex plane must be known to calculate the quark distribution. The present investigation enlarges the class 
of model propagators which can be used to calculate light-cone dominated amplitudes. We show how meromorphic 
propagators can be used in Minkowski space to arrive at parton and generalized parton distributions as well as the 
electromagnetic form factor expressed on the light cone. Additionally vertex functions for space-like processes can be 
modeled using meromorphic functions which allows one to go beyond the impulse approximation. 

Complex conjugate singularities present in solutions to Dyson-Schwinger equations have been studied in the con- 
nection with the violation of Osterwalder-Schrader reflection positivity and confinement || . Recent work [|J in solving 
the Bethe-Salpeter equation with a quark propagator consisting of pairs of complex conjugate singularities shows that 
the width for meson decay (into free quarks) generated from one pole is exactly canceled by the contribution from its 
complex conjugate. Additionally recent studies have modeled Euclidean space lattice data with propagators that have 
time-like complex conjugate singularities [13, E|- I n this work, we pursue the calculation of space- like amplitudes in 
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Minkowski space using a simple model propagator consisting of a pair of complex conjugate poles. Specifically, we 
are interested in the calculation of light-cone dominated amplitudes for this model which necessitates a treatment in 
Minkowski space. 

The paper is organized as follows. First in section [n] we present the scalar model used to investigate complex 
conjugate poles in propagators. Here we show that naive calculation of the quark distribution in Minkowski space 
is problematic. The distribution has neither proper support nor is positive definite. In section ITTT1 we demonstrate 
that despite troubles in Minkowski space, amplitudes involving propagators with pairs of complex conjugate poles 
can be calculated directly in Euclidean space. The double distribution (DD) is extracted from non-diagonal matrix 
elements of twist-two operators. This distribution satisfies all of the relevant spectral and support properties and 
can be used to derive parton and generalized parton distributions for the model, in addition to the electromagnetic 
form factor. In section IIVI we show that the quark distribution can be calculated in Minkowski space by using a 
natural modification of the cutting rules applied to the handbag diagram. The effect of these cutting rules is to put 
intermediate states on their complex mass shells. Finally in section we rectify the situation in Minkowski space 
by analytically continuing amplitudes from Euclidean space. This justifies the cutting rules presented. The detailed 
calculation of the generalized parton distribution in Minkowski space is contained in the Appendix. After analytic 
continuation, the model's parton and generalized parton distributions agree with those calculated from the Euclidean 
space DD. A brief summary (section lyijl concludes the paper. 



II. PROBLEMS IN MINKOWSKI SPACE 



In this section we present the simple model under consideration. Here the propagator is treated in Minkowski space 
where difficulties are encountered. We show that amplitudes cannot be directly calculated in Minkowski space when 
complex conjugate poles are present in propagators and vertices. 

The model we take is </> 3 theory with electromagnetic interactions. Equivalently we can view this model as a bound 
state of two scalar particles with a trivial Bethe-Salpeter vertex T(k, P) = 1, where the coupling constant is assumed to 
be absorbed into the overall normalization. We make a simple Ansatz for the non-perturbative propagator consisting 
of a pair of complex conjugate poles 

s(k) = ^ 2 -a 2 + b 2 ) m 

{ ' {k 2 -a 2 + b 2 ) 2 + Aa 2 b 2 ' 1 ' 

where a 2 — b 2 > 0. Defining for ease m 2 = a 2 — b 2 and e = 2ab (which is taken to be positive without any loss of 
generality), we can write the propagator as 

<SY/c) = — — ■ - ; -)- — ~ ■ - ; . (2) 

K — m z + ie k z — m z — it 
The light-cone energy poles 1 of the propagator are thus 

_ k^ 2 +m 2 _ ie^ 
a ~ 2k+ 2k+ [6) 

and k~, = where * denotes the complex conjugate. Although we use a scalar model, results straightforwardly 

extend to spin-^ particles, e.g., since only the pole structure of Eq. Q is relevant. 

Now let us consider calculating this model's quark distribution by projecting onto the light cone. The quark 
distribution can be derived by fixing the plus- momentum of the active quark x = k + / P + , see Figure ^ and taking 
the plus-component of the current. Thus up to overall normalization, we have the expression 

q{x) cx J d*k 5{k + - xP + ) xS{k)S{k - P)S(k). (4) 

The k~ integral is then performed by residues. Choosing a frame in which = 0, the spectator propagator has 
light-cone energy poles 

. I 2 o 

6 + 2(k+ - p+) ~ 2(k+-p+y (5) 



1 For any vector a M , we define the light-cone variables = -h= (a ± a 3 ) . 



3 




FIG. 1: Triangle diagram at zero momentum transfer used to calculate the quark distribution function by projecting onto the 
light cone. 



and k b , = (k b )*. 

Performing the fc~ integral in Eq. Q, we arrive at the quark distribution 

) [ Res(fc~, ) + Res(fc b ". )] + 6[x(l - x)} [Res(fc~.) + Res(fcj^)] + 6{x - 1) [Res(fc~„) +Res(fc fe ;)]) (6) 

This distribution does not have proper support, i.e. it is non- vanishing outside the interval x € [0, 1]. Moreover, the 
distribution is not real valued, whereas it should be positive definite. Thus the model based on the propagator in 
Eq. m cannot be suitably formulated in Minkowski space. We will find below that Eq. (JTJ makes sense as a Minkowski 
space propagator only after analytic continuation from Euclidean space for the amplitude in question. 



III. COVARIANT CALCULATION IN EUCLIDEAN SPACE 



In Euclidean space, the model propagator is 

S E (k) = Y ,„ 1/2 9 — -. (7) 

£=± 

Here and below we use the shorthand e = ± to denote the pair of poles e = —2ab, +2ab. Unlike in Minkowski space 
where the measure is imaginary, contributions to Euclidean space amplitudes are real and one has no difficulty in 
calculating form factors and distribution functions using Se (k) in the relevant diagrams. The simplicity of the model 
at hand will allow us to calculate its double distribution analytically and thereby determine the quark distribution 
and electromagnetic form factor, since these functions are related to the double distribution by the so-called reduction 
relations. The remainder of the paper will be devoted to calculation of these quantities in Minkowski space by 
projecting onto the light-cone. 

GPDs are not Lorentz invariant objects, however, they stem from a projection of a Lorentz invariant double 
distribution function ^4|. These functions are particularly attractive from the perspective of model building [l5| . 
though one must be careful that the starting point is indeed covariant , otherwise desirable distribution properties 
and straightforward physical interpretation may be sacrificed. The model under consideration is fully covariant, and 
thus the DD representation is an ideal testing ground for our model propagator. Hence we proceed to calculate the 
model's Euclidean space DD, recalling along the way the relevant properties of DDs. 

Let = <9 M — d^. For this scalar model, we define the twist-two operator of spin-n as 

O w -''» = <)>(0)iD^iD^ ■■■W^i (f>(0), (8) 

where the action of I ■ ■ • 1 on Lorentz indices produces only the symmetric traceless part. 

We work in Radyushkin's asymmetric frame 2 with P as the momentum of the initial state, P + A that of the final 
and t = A 2 . Following the two-component formalism of the non-diagonal matrix element of O m '"' i " can be 



2 The non-diagonal matrix elements of twist-two operators are, however, more conveniently expressed in variables symmetric with respect 
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k + A 



P P + A 

FIG. 2: Diagram used to calculate non-diagonal matrix elements of twist-two operators (denoted by a cross) 



decomposed into Lorentz invariant moment functions A nk (t) and B nk {t) 

r " i 

(P + A\0^-^\P) = (2P + A)^V — - ,.. A nk (t) (2P + Af 1 •••(2P + A)^"- fc (-A) 

z — ' kl(n — k)\ 



k=0 



fln-k + l . . . 



n 
fc=0 



fc!(n-fc)! 



P„ fc (i) (2P + A) P1 • • • (2P + A)^- fc (- A) 



(9) 



Hcrmiticity forces the matrix elements of O fJ,,Ml "' IJ ' n to be invariant under the transformation 

jp -> P + A 
|A^-A 

Consequently the values of are restricted to be even in the first sum and odd in the second. As it stands there 
is considerable freedom in this decomposition, e.g. one could rewrite the above with kB n ^-i(t)/(n — k + 1) as a 
contribution to A nk {t). Carrying this out for all k, puts the bulk in the first term and renders the second term 
proportional only to the symmetric traceless part of (n + 1) A's — moments of the Polyakov- Weiss £>-term 19] . This 
is the usually encountered form of the DD with D-teim. Alternatively one can also express the moments as projections 
of a single Lorentz invariant function [2(J. Calculationally, however, Eq. @ becomes the most practical to work with 

na. 

The F and G DDs can be defined as generators of the coefficient functions 



A nk {t) 
B n k{t) 



1 rl-X 

dx / ^ 
i) Jo 

1 pl-x 

dx / 
i) Jo 



dy x n ~ k (x + 2y~l) k F(x,y;t) 
dy x n ~ k (x + 2y-l) k G(x,y;t). 



(10) 
(11) 



As a consequence of the restriction on k in the sums, the function P(x, y; t) is Miinchen symmetric |'2'2t . i.e. F(x, y; t) — 
F(x, l — x~y;t), while G(x, y; t) is Miinchen antisymmetric. Also for n-even, there is no contribution from the D-term 
to the function G{x,y;t). 

These functions then appear in the decomposition of matrix elements of light-like separated operators 



( P + A | 0(0) iz-D <p(z~) | P ) 

"1 rl — X pi pl — X 

dx I dy er txP - z+lyA - z F{x, y; t) - A ■ z / dx dy e -«^+n/A- z 



(2P-Z + A-2) 



II 



II 



II 



G(x,y;t), (12) 



o 



to initial and final states. Since we use perturbative diagrams and do not have antiparticles, asymmetrical variables are warranted. 
Good discussion of the conversion from symmetrical and asymmetrical variables and distributions can be found in Il7l . Additionally 
advantages and disadvantages of both are presented. 
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where z 2 = 0. 

Denoting Q = — A+/P+ > 0, the GPD in asymmetric variables reads 

r J _- ixP + z~ <-> 

ff(x,c,*) = y 27r(2 _ c) (p+&\m*D + <f>(z-)\p)- (13) 

Physically C plays the role of Bjorken variable for deeply virtual Compton scattering. Inserting Eq. (|12f) into this 
definition yields 

H(x,(,t)= [ dz f dyS(x-z-Cy)\F(z,y;t) + -^-G(z,y;t)]. (14) 
Jo Jo 1 1 s J 

By integrating Eq. (|14f> over x, we uncover two sum rules: the sum rule for the form factor 

/ dx / dy F(x,y;t) = F(t) (15) 
Jo Jo 

and the G-sum rule 

dx dyG(x,y;t) = 0, (16) 



which follows since G is Miinchen antisymmetric. Eq. Ijl6|l is important and mandated by current conservation. All 
too frequently the G DD function is overlooked and treated as identically zero. Lastly, the quark distribution function 
q{x) can be found from the DD at zero momentum transfer, see Eq. (|1H|) . 

q(x) = / dyF(x,y;0). (17) 
Jo 

We can use the decomposition in Eq. to calculate our simple model's DD. Parameterizing the momenta as in 
Figure |3 the non-diagonal matrix element of 0^ n > reads 

(P+ MO»"»*n\P) = ™ V fd% (2fc + A)^(2fc + A)^--(2fc + AH 

1 17 tt2 / [fea + m 2-i e ][(fc + A) a + m a -i^[(fe-P) a +m a -ie / '] 



7T , „ 

e,e',e"=± 



The normalization constant iV is chosen by the condition F(0) — 1. Let us denote the propagators simply by 
21 = (k — P) 2 + m 2 — i e", 03 = (k + A) 2 + rn 2 — i e' and £ — k 2 + m 2 — it. We introduce two Feynman parameters 
{x,y} to render the denominator specifically in the form [a;2l + j/23 + (1 — x — y)£]~ 3 . One then translates k^ to render 
the integral (hyper-) spherically symmetric via the definition k^ = + xP^ — yA. 11 . The resulting integral over / can 
be evaluated directly (remember we are in Euclidean space). 

Binomially expanding the result of the integral, we can make contact with Eq. JUJ) and subsequently determine the 
F and G double distributions by inspection from Eqs. (|10fl and l|ll|) . Defining the auxiliary functions 

D (x, y; t) = m 2 - x(l - x)M 2 -y{l-x- y)t (19) 

and 



D (x,y;t) 

'\x,y-,i*) = -Li ? —, — 

r. P>o(X,1 

z—0/x,y,x-\-y v L 

the DDs can be written simply as 



D(x,y;t) = N ]T T2 9V » ' ( 20 ) 



F(x,y;t)=xD(x,y;t) (21) 
G(x,y;t) = (x + 2y-l)D(x,y;t). (22) 

Accordingly F is Miinchen symmetric and G is antisymmetric. Notice although e is finite, corresponding results using 
the standard perturbative propagator can always be recovered in the limit e — > 0. For example, the correct F and G 
DDs are recovered in the limit e — » l21l. 
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FIG. 3: On the left, the quark distribution Eq. is plotted as a function of x for a few values of e in GeV 2 . On the right, 
the form factor calculated from Eqs. 1151 and I2H is plotted as a function of —t for a few values of e in GeV 2 . The model 
parameters are arbitrarily chosen as: M = 0.14 GeV and m = 0.33 GeV. 




FIG. 4: Plots of the GPD calculated from Eqs. ||3T),||32J and itTHl . On the left, the GPD is plotted as a function of x for a few 
values of e (in Gev 2 ) at fixed £ = 0.5 and t = —1.0 GeV 2 . On the right, the GPD appears at fixed e = 0.17 GeV 2 and t = —2.0 
GeV 2 and is plotted as a function of x for a few values of f. The model parameters are arbitrarily chosen as: M = 0.14 GeV 
and m = 0.33 GeV. 



The model GPD can be derived by utilizing Eq. 1|14JI . although the integral must be performed numerically. The 
quark distribution can be found via the reduction relation Eq. (|17(l . namely 

a(x) - N V ( *(1-*)A,(*.0;0) x t e(l-2z) \ 

~ \ D ° {x > 0; 0)2 +e2{1 - 2z)2 e D ° {x > 0; 0) / 

Lastly the form factor can be found from the sum rule Eq. (|15|) 

In Figure 3, we plot the quark distribution and electromagnetic form factor for various values of e in GeV 2 . We 
have arbitrarily chosen the other model parameters as M = 0.14 GeV and m — 0.33 GeV. Additionally in Figure 4, 
the GPD is plotted: first at fixed £ and t for various values of e and then at fixed t and e for various values of £. 
Curves corresponding to e = are the standard results for a propagator with one real pole. 



IV. CUTTING RULES 



In this section, we show there is some hope in working with the model propagator Eq. JQ| in Minkowski space. We 
demonstrate that the quark distribution can be derived by a straightforward generalization of the cutting rules. 

Consider the forward Compton amplitude Ai^ depicted in Figure |SJ In the Bjorken limit, the imaginary part 
of this diagram is related to the quark distribution. For simplicity, we can choose a frame in which = 0. The 
minus-plus component of the forward Compton amplitude in such frames reads 

iM~ + = J d*k {2k- + q-)S(k) S{k - P) S{k) (2k+ + q + )S{k + q). (24) 

In the scalar particle case, the minus-plus component of the forward Compton amplitude can be used to define 
the quark distribution in a way analogous to the spin-i case. The relation is simply 3(A4 ^) cx q(x), see |23j . 
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P P 

FIG. 5: Handbag diagram for the forward Compton amplitude. Dashed line denotes the cut which yields the quark distribution 
in the Bjorken limit. 



In Eq. ()24fl. we have adjusted the overall normalization so that equality between 3(.M~ + ) and q(x) holds. With 
standard perturbative propagators, we could calculate the imaginary part of M.^ v by using the cutting rules, whereby 
one replaces the cut propagators in Figure [5] by an on-shell prescription, namely 

S{k) e=0 — > -2mS(k 2 -to 2 ). (25) 

In Eq. (|25[l . we have specified e = for the case of the free particle propagator. Since large momentum flows through 
the handle of the handbag, we may neglect the mass of the struck quark and use the standard cutting rule for S(k + q). 
In the Bjorken limit, we define x = —q 2 /2P ■ q which remains finite as q 2 , P ■ q — > oo and is kinematically bounded 
between zero and one. Further we orient our frame of reference so that q has a large minus component in this limit, 
and consequently q + = q 2 /2q~ is finite. Hence we have the familiar replacement 

S(k + q)— -> 5{k+-xP+). (26) 

2q- 

To complete the cut, we must deal with the spectator particle's complex mass shells. We must worry about the 
propagator Eq. Q only where the denominator is zero. Thus we are lead to the cutting rule for the propagator Eq. Q 



S(k) 



S(k 2 - m 2 + i e) + S(k 2 



ie) 



(27) 



which puts the intermediate state on its complex mass shells. Furthermore, the limit e — > produces the regular 
cutting rule Eq. 

Using this cutting rule for the spectator particle along with Eq. I|26|). we can deduce the quark distribution from 
Ss(Ai h ) in the Bjorken limit 



q{x) 



AN 

7T 



d 4 k 5(k+ - xP + ) S(k- - fcr) + 5(k~ - kZ) 



xS(k) 2 
l-x 



(28) 



where the light-cone energy poles are given in Eq. JSJ. Notice the resulting distribution is real and has proper support 
due to the kinematic constraint x G [0, 1] imposed by the Bjorken limit. Evaluation of the two trivial integrals leaves 
only the transverse momentum integration 



q(x) = 



N 



E 



dk 1 



x(l — x) 



B w (x, k- 1 , e, e" \M 2 ) D w (x, , e' , e" \M 2 ) 



(29) 



where we have defined the Weinberg propagator generalized for complex masses as 



D w (a;,k- L ,e,e'|M 2 



M 



x(l — x) x 1 — x 



(30) 



Evaluation of the k 1 - integral yields an analytic expression for q(x), which is identical to that obtained from the DD 
Eq. H23|) . Notice Eq. (|28|l is equivalent to evaluating Eq. at Kes(k^) + Res(k^,) and hence the quark distribution 
is derived as if the complex conjugate spectator poles both lie in the upper-half complex plane! We will understand 
this better once we analytically continue from Euclidean space. 
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x < 



M 



x > 



Re(kT ) 



Re(k- ) 



FIG. 6: Complex light-cone energy plane for the propagator. 



V. ANALYTIC CONTINUATION 



Above we have seen that modified cutting rules can be used to derive the correct quark distribution function in 
Minkowski space. In essence the result stems from putting the spectator particle on its complex mass shells. This will 
be justified by careful analytic continuation of Euclidean space amplitudes. Below we consider the Minkowski space 
calculation of the generalized parton distribution which cannot be derived by cuts. This will force us to deal with the 
underlying Wick rotation necessary to define the model in Minkowski space. 

In section IITT1 the model double distributions were calculated in Euclidean space. Thus to calculate related ampli- 
tudes in Minkowski space, we must Wick rotate in the complex energy plane: k 4 — > ik°. The analytic continuation 
can be viewed alternately in the complex light-cone energy plane. For k 3 = 0, the rotation is from the axis to 

the 3?(fc _ ) axis. In the general case, such a correspondence can only be made precise by considering the Wick rotation 
in terms of k° and then boosting to the infinite momentum frame. This requires tedious algebraic manipulations and 
a proliferation of energy poles and time-ordered diagrams. Indeed it is easier just to imagine the rotation simply and 
deal with the light-cone singularities. This is the approach we present. 

Before tackling generalized parton distributions in Minkowski space, let us imagine a simpler fictitious example. 
Consider some well-defined Euclidean space amplitude having only poles at k~ and k~, , see Eq. Q . To calculate the 
amplitude in Minkowski space, we naively integrate along the ) axis. In general the correct path on which to 
integrate is one which nears the ^R(kr) axis except for detours around energy poles in the first and third quadrants. 
Such a path is correct since it can be continuously deformed into the Euclidean path. The difference between the naive 
integration and the correct path is a sum of residues of the Wick poles. The energy poles of our fictitious amplitude 
are depicted in Figure HO Their location depends upon the sign of x = k + / P + . Thus for this amplitude the correct 
continuation from Euclidean space is 

J dSf(k~) — ► j (Mt{k-) + 2m[6(-x) - 9(x)] Res(fc^)- (31) 

Notice only k~* is a Wick pole; this is expected because we know the limit e — > can be analytically continued in 
the naive fashion. Closing the contour in the upper-half plane to perform the Minkowski space energy integral 3 and 
evaluate our fictitious amplitude, we pick up 2m[6{— x) Res(k~) + 9(x) Res(fc~„ )]. The net result according to Eq. (|3T|) 
is thus 

2m6{-x) [Res(fc-) + Res(fc~*)] • 

Looking back at Figure El the net result after Wick rotation amounts to both poles lying in the same half-plane; or 
equivalently, we have effectively integrated in either the right- or left-half plane. 



3 One must be careful of zero modes [248 for which fc + = 0. In such cases, the pole lies on the contour at infinity and the integration 
cannot be performed by residues. Since our fictitious example is only schematic, we are neglecting the issue of zero modes and are hence 



excluding amplitudes of the form 



d 4 k 



which must be handled separately. The example 



(k 2 - m 2 +ie) T ' 
d 4 k 



n > 2, 



(k 2 — m 2 + i e) 2 (k 2 — m 2 — i e) 2 



is devoid of zero-mode complications and more closely parallels the expressions encountered for GPDs. 
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The space-like amplitude 4 for the generalized parton distribution can now be continued to Minkowski space and 
hence be evaluated by projecting onto the light cone. To do so, we refer to Figure [21 and insert the non-local light- 
cone operator <f)(0)iD + 4>(z~) in place of the local twist-two operators denoted by a cross in the Figure. Here the 
plus-component picks out the leading-twist contribution according to light-cone power counting. Thus in momentum 
space, we arrive at 

H(x, C, t) = 2N lf ln I d 4 k 8(k+ - xP+) (2k+ + A+) S(k)S(k - P)S(k + A). (32) 



l-C/2 

Above we have included the ^-dependent pre-factor to normalize the action of D + between non-diagonal states. The 
overall normalization is then the same as in Eq. I|18fl . By writing Eq. i|32|) in Minkowski space, we must also keep in 
mind the Wick residues implicitly necessary so that Eq. (|32[l is meaningful. In addition to the poles k~ , k^ (given in 
Eqs. © and J5J), respectively) and their complex conjugates, the integrand of Eq. (|32[1 also has the poles 

(k ± +A ± ) 2 +m 2 it 
c " " + 2(fc+ + A+) 2(fc+ + A+) (33) 

and k~* = (k~)*. Carrying out the light-cone energy integration in Eq. I|32(l as well as adding relevant residues 
resulting from the Wick rotation produces (the subtle details of this calculation appear in the appendix) 

H (x, C, t) = -2m 9(x)6(C - x) \ Res(fc") + Res(fc Qt )1 + 2m 9{x - £)6{l - x ) \ Res(fc^ ) + Res(k^ ) , (34) 

where the residue is of the integrand in Eq. (|32|l . As a result of the effective relocation of poles to the same half-plane 
as their complex conjugates, the resulting GPD Eq. (|34|l is real and vanishes outside x from zero to one. 

Using the Weinberg propagator Eq. Ij30|l . the residues can be compactly written in terms of relative momenta. 
Defining the relative momentum of the final state as 



y _ x ~ C 

" i-C 

and the relative momentum of the photon as 



(l-x')A^, (35) 



x 

c 



k" 1 - = k ± + x" A ± , (36) 

the light-cone GPD can be expressed in the form 

(1 - C/2) H(x, C, t) = 8(x)6(( - x) H, (x, C, t) + 6(x - 06(1 - x) H 2 (x, C, t), (37) 
where we have made the abbreviations 

H 1 (x,C,t) = (2*"-l)£- J x „ (1 _ d X)(l - x) Dw{x ' k± > £ ' £ " |M2) Dw(x "' k,,± > £ ' £ ' 1 l) (38) 

e,e',e" — ± 

H 2 (x,(,t) = (2x-0^ E /^^w D w(^k ± ,e, e ''|M 2 )Dw(x',k' ± ,6',e''|M 2 ). (39) 

e,e',e" — ± 

Firstly one can see analytically that the correct quark distribution results at zero momentum transfer, namely 
H 2 (x,0, 0) = q{x), where q(x) is given by Eq. lt2P|) . As remarked in section llVl this function is identical to that 
obtained from the double distribution Eq. iffifl) . Secondly, the resulting light-cone GPD Eq. l|3T|l agrees numerically 
with that found from the double distribution, via Eq. I|14f) . which is plotted in Figure 4. Finally the electromagnetic 
form factor found from the sum rule 



F(t) = f 
Jo 



dxH(x,C,t) (40) 



4 There are additional complications for time-like amplitudes and for amplitudes involving unstable bound states. In these cases Wick 
poles are present even when standard perturbative propagators are used. The analysis above must be more carefully considered in these 
cases where threshold effects arc already inherent in the analytic continuation to, or from, Euclidean space. 
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agrees numerically with the result of Eq. (|15|l (which is plotted in Figure 3). The £-independence of Eq. I|40|) stems 
from Lorentz invariance which is present, however, not manifest in Eq. i|37|) . Thus with the calculation of Eq. i|37|) 
from analytically continuing to Minkowski space, space-like amplitudes now agree with those calculated from the 
Euclidean space double distribution. 

VI. SUMMARY 

Above we consider calculation of amplitudes for space-like processes using a scalar propagator with one pair of 
complex conjugate poles. Although we use a scalar model, generalization to higher spins is clear since the energy 
denominators are universal. Moreover the analysis can be extended easily to the case where vertex functions have com- 
plex conjugate singularities. Such models cannot be directly employed in Minkowski space, they must be analytically 
continued from Euclidean space. 

In section [n] the problems of using a propagator with complex conjugate poles in Minkowski space are discussed 
at the level of the quark distribution function. If the model is defined in Minkowski space, one will generally violate 
the support and positivity properties of the quark distribution. Next in section IIIII we show the model is perfectly 
well defined in Euclidean space by calculating non-diagonal matrix elements of twist-two operators. This leads us to 
the model's double distribution which we used to calculate partem and generalized parton distributions as well as the 
electromagnetic form factor. 

In the remainder of the paper, we investigate how to calculate amplitudes properly in Minkowski space. Firstly 
amplitudes dependent on the imaginary part of some set of diagrams can be calculated by using a straightforward 
generalization of the cutting rules. We apply this to calculate the quark distribution function from the handbag 
diagram in the Bjorken limit in section llVl Lastly we consider the analytic continuation of space-like amplitudes from 
Euclidean space. This is complicated by the presence of Wick poles and requires their residues to be appropriately 
added when amplitudes are calculated. The details of the GPD calculation appear in the Appendix. Resulting 
functions calculated in Minkowski space after the Wick rotation agree with those obtained from the model defined in 
Euclidean space. 

This work lays the foundation for calculating light-cone dominated amplitudes using meromorphic model propaga- 
tors and vertices constrained by lattice data 5 and Ward-Takahashi identities. Distribution functions for such models 
could be calculated rigorously since the pole structure of the propagator is known and relevant integrals converge in 
the complex light-cone energy plane. As far as light-cone phenomenology is concerned, resulting expressions would 
be truly Poincare covariant (as opposed to diagonal with respect to the non-interacting operators) and would satisfy 
field-theoretic identities. Filling these two gaps is essential for adequate hadronic phenomenology for processes at 
large momentum transfer. Moreover, such models would help light-cone methods and Dyson-Schwinger studies reach 
complementary standing. 
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APPENDIX: CALCULATION OF THE GPD 

Below we derive Eq. (|34() for the GPD in Minkowski space. The details have been relegated here since there is 
some subtlety. In order to evaluate Eq. l|32fl . which implicitly needs analytic continuation, we must shift the energy 
integration variable and define a prescription for dealing with vanishing real parts. Let us see how these difficulties 
arise. 

In considering the Wick rotation, one is usually only concerned with the single denominator that results from 
combining propagators via Feynman parameters. For the moment, let us ignore the complication of complex conjugate 
pairs of poles. In this case, combining the denominators of Eq. I|18f) using Feynman parameters results in [I 2 — 



One must proceed with caution: the lattice calculations employ Landau gauge, while we have tacitly used light-cone gauge above. 
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FIG. 7: Complex light-cone energy plane for the shifted poles Eq. 1A.3I for the generalized parton distribution. 



D (x,y;t) + «e]~ 3 , where D is given by Eq. I|19(l . Since the bound state is stable and t is space-like, D (x,y;t) is 
always positive and hence there are no Wick poles. Moreover, we have shifted the variable k^ to arrive at l* 1 . 

In analytically continuing the expression with uncombined propagators (and again no complex conjugate poles), 
we are confronted with a problem. The pole k^ , for example, is shifted by the energy P~ . Thus the location of the 
singularity in the complex plane will be shifted parallel to the real axis depending upon the relative magnitude of the 
spectator's kinetic energy and P~ . In the schematic example shown in Section the pole k~ does not have such a 
shift. Thus for the b and c poles, the location of the singularities depicted in Figure will shift along the real axis 
(depending on P~ and A - ) and there will be threshold Wick poles [the threshold is defined when ^ft{k b ' c = 0]. This 
is unphysical: we just demonstrated the Wick rotation can be done for the combined denominators without crossing 
any poles. To perform the same Wick rotation at the level of separate propagators, we must use the freedom to shift 
the energy variable as well as the stability of the bound state. 

On the light-cone, the bound state stability condition can be expressed as 



where k re i = — xP 1 - is the relative transverse momentum of the two constituents. Because we consider the elastic 
electromagnetic form factor, there is an analogous relation for the final state P 1 . Since each propagator contains the 
kinetic energy of a single particle, the bound state stability condition can never be utilized without shifting k~ . Yet 
in order to perform such a shift, the real part of one pole must be zero and hence we must invent a prescription for 
moving this pole off the Euclidean contour. 

We choose the translation: k~ — > k~ + 5k~ + rj, where 77 is a positive infinitesimal and 

The resulting poles of the integrand in Eq. (|32|l we denote 

1 = K - <> k = - P + 2P+i{l-a;) ~ IP^h 

K = K -Sk--r)= -r] - 2P+ }^-i) • ( A -3) 

~K = K - 5k- = -P'~ 



2P'+x'(l-x') 2P+(x-Q 

and similarly for their complex conjugate partners. Notice imaginary parts of the poles are unaffected by the energy 
translation. The r\ prescription has displaced the resulting spectator pole away from the Euclidean path independent 
of x. Moreover the a and c poles have non-zero real parts; so ry has been set to zero for these poles above. 

Using the expressions for the new poles Eq. I|A.3|) and the bound state stability condition Eq. I|A.1|) . we can 
determine the quadrant location of the singularities independent of P~ and A - (see footnote 4). These quadrant 
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locations are depicted for the full range of x in Figure Accordingly poles in the first and third quadrants are Wick 
poles. As we saw in Section the net result of the analytic continuation is to evaluate the integral by effectively 
closing the contour in the right- or left-half plane. Hence the GPD Eq. (|32|) is 



H(x,C,t) = -2irie(x)9(C-x) Res(fc~) + Res(fc"«) + 2m 8(x - £)0(l - x) Res(fc") + Res(fc") 



(A.4) 



Evaluating the residues in Eq. (|A.4() yields the result of Section [VJ namely Eq. I|37|) which is algebraically equivalent 
to Eq. PI). 

Notice from Figured other t\ prescriptions for the shift, such as + 2 p+(x-i) ' ^ eac ^ ^° an mcorrec t result for the case 
when there are no complex conjugate pairs. The figure shows that the infinitesimal prescription must be positive and 
independent of the sign of a;, x — £, etc, in order to reproduce the familiar result. It is interesting to note that for 
x > 1 in the case without conjugate pairs, all poles of the integrand are Wick poles. Though since the integral is 
convergent, the sum of these Wick residues vanishes. 

The interested reader can verify that the alternate shifts which use the same pole prescription 



2P+3 



V 



(A.5) 



2P+(x-() 



also yield the correct results provided t is space-like. 
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